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Abstract — A  compact  model  and  simulation  methodology  for 
chalcogenide  based  memristor  devices  is  proposed.  From  a 
microprocessor  design  view  point,  it  is  important  to  be  able  to 
simulate  large  numbers  of  devices  within  the  integrated  circuit 
architecture  in  order  to  speed  up  reliably  the  development 
process.  Ideally,  device  models  would  accurately  describe  the 
characteristic  device  behavior  and  would  be  represented  by 
single-valued  equations  without  requiring  the  need  for  recursive 
or  numerically  intensive  solutions.  With  this  in  mind,  we  have 
developed  an  empirical  chalcogenide  compact  memristor  model 
that  accurately  describes  all  regions  of  operations  of  memristor 
devices  employing  single-valued  equations. 

Keywords-  Memristor;  Memory  Device;  Thin  Film  Device 

I.  Introduction 

The  memristor  device  postulated  in  1971  by  Leon  Chua  [1] 
as  the  fourth  basic  circuit  element  has  received  much  attention 
in  the  research  community  since  the  publication  of  Strukov’s 
2008  paper  titled  “The  missing  memristor  found”  [2]  The 
memristor  name  is  a  contraction  for  memory  resistor  [1] 
because  that  is  exactly  its  function:  to  remember  its  history  [3]. 
The  memristor  is  a  two  terminal  passive  device  whose 
resistance  state  depends  on  its  previous  state  and  present 
electrical  biasing  conditions,  and  when  combined  with 
transistors  in  a  hybrid  chip,  memristors  could  radically  improve 
the  performance  of  digital  circuits  without  the  need  for  farther 
reduction  of  transistor  dimensions  [3].  Given  their  two  terminal 
structural  simplicity  and  electronic  passivity,  the  applications 
for  memristor  technology  range  from  non-volatile  memory, 
instant  on  computers,  reconfigurable  electronics  and 
neuromorphic  computing  [3,4].  According  to  Chua  [4],  the 
memristor  behaves  like  a  linear  resistor  with  memory  but  also 
exhibits  many  interesting  nonlinear  characteristics,  and  several 
electronic  models  have  been  presented  to  describe  the  electrical 
behavior  of  memristor  devices  [1,2, 4-6].  However,  given  that 
memristor  devices  are  not  commercially  available,  good 
physical  model-to-hardware  correlations  have  not  been  yet 
been  reported  in  the  published  literature.  Therefore,  in  this 
work,  we  present  what  we  believe  to  be  the  first  model-to- 
hardware  correlation  of  memristor  electrical  characteristics.  In 
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our  studies,  we  have  employed  both  linear  and  nonlinear 
models  from  the  published  literature  to  fit  our  memristor 
hardware.  However,  we  have  observed  that  these  published 
models  do  not  represent  accurately  the  electrical  characteristic 
behavior  or  our  memristor  device  hardware.  Therefore,  we 
have  developed  a  simple  compact  model  that  accurately 
represents  the  electrical  behavior  of  chalcogenide  based 
memristors.. 

II.  Experimental  Details 

Memristor  devices  were  fabricated  on  200  mm  p-type  Si 
wafers.  Isolated  tungsten  bottom  electrodes  were  patterned  on 
the  wafers  and  a  planarized  nitride  layer  was  used  for  device 
isolation.  Vias  were  etched  through  the  nitride  layer  to  provide 
contact  to  the  bottom  electrode  and  to  define  the  device  active 
region.  Prior  to  deposition  of  the  memristor  materials,  the 
wafers  received  an  Ar+  sputter  clean  to  remove  residual 
material  and  potentially  remove  any  oxides  that  might  have 
formed  over  the  tungsten  electrode.  The  memristor  device 
structure  consists  of  the  layers  (from  bottom  electrode  contact 
side  to  top  electrode  contact):  300  A  Ge2Se3/500  A 
Ag2Se/100  A  Ge2Se3/500  A  Ag/100  A  Ge2Se3.  The  100  A 
Ge2Se3  layers  are  needed  for  device  processing  only  since  Ag 
cannot  be  deposited  directly  on  Ag2Se  and  since  tungsten  (for 
the  top  electrode)  does  not  adhere  well  to  Ag  in  this  material 
stack.  The  layers  were  deposited  by  thermal  evaporation  using 
a  CHA  Industries  SE-600-RAP  thermal  evaporator  equipped 
with  three  200  mm  wafer  planetary  rotation.  The  rate  of 
material  deposition  was  monitored  using  an  Infineon  IC  6000 
with  a  single  crystal  sensor  head.  The  base  system  pressure  was 
1x10-7  Torr  prior  to  evaporation.  A  tungsten  top  electrode  was 
deposited  by  sputtering  tungsten  (350  A)  and  etching  to  define 
the  device  top  contact  and  bond  pad.  Etching  was  performed 
with  a  Veeco  ion-mill  by  etching  through  the  tungsten  and  the 
memristor  device  materials  and  stopping  on  nitride.  The  top 
and  bottom  electrode  bond  pad  contacts  were  80  um  x  80  um. 

Electrical  measurements  were  performed  with  either  an 
HP4145B  or  an  Agilent  B1500A  semiconductor  parameter 
analyzer.  A  Micromanipulator  6200  microprobe  station 
equipped  with  a  temperature  controllable  wafer  chuck  was  used 
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for  the  wafer-level  device  measurements.  Probe  tips  were 
either  Micromanipulator  W  size  7B  or  American  Probe  and 
Technologies  73CT-APTA  probe  tips.  The  tested  devices  were 
180  nm  in  diameter  as  defined  by  the  via  etched  through  the 
nitride  layer  to  the  bottom  electrode.  As  fabricated,  the  devices 
do  not  display  negative  differential  resistance  (NDR).  NDR  is 
induced  in  the  devices  by  application  of  short  (ns)  negative 
potential  pulses  (negative  with  respect  to  the  top  electrode). 
Application  of  several  of  these  pulses  irreversibly  induces  NDR 
in  the  devices  allowing  them  to  permanently  function  as 
memristor  devices. 
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Fig.  2.  Memristor  device  modeling  structural  diagram. 

III.  Memristor  Models 

The  most  basic  mathematical  definition  of  a  memristor  is 
that  of  a  current-controlled  device  for  circuit  analysis  in  the 
generalized  class  of  nonlinear  dynamical  systems  called 
memristive  systems  described  by  the  equations  [2,4] 


«  »  i  i  i  i  •  i  i 


•0.5  -0.4  -0.3  -0.2  -0.1  0  0.1  0.2  0.3  0.4  0.5 

Input  Voltage  (V) 


Fig.  1 .  (a)  Memristor  fabrication  design  and  (b)  typical  memristor  device 
Lissajous  I-V  curve  characteristic  behavior  atlOOHz  frequency. 

Figure  1(a)  displays  the  memristor  typical  device  structure, 
and  figure  1(b)  shows  the  electrical  Lissajous  I-V  curve 
characteristic  response  behavior  of  a  memristor  device  under  a 
sinusoidal  input  of  0.5  V  volts  amplitude  and  100  Hz 
frequency.  From  the  figure  we  can  clearly  observe  that  the 
memristor  device  toggles  between  two  states  of  high  and  low 
conductivity.  As  the  device  transitions  between  low  to  high 
conductivity  states  it  goes  through  high  nonlinear  diode-like 
processes  at  approximately  -0.35  and  0.2  V  threshold  voltages 
respectively.  The  threshold  voltage  analogy  is  used  here  to 
describe  the  biasing  regions  where  nonlinear  behavior  occurs. 


v  =  R(w,  Qi 


(1) 


dw 

dt 


=  f(w,  0 


(2) 


where  w  can  be  a  set  of  state  variables  and  R  and  f  can  in 
general  be  explicit  functions  of  time  [2],  [4].  For  simplicity  and 
ease  of  simulation,  the  memristor’ s  resistance  or  memristance 
definition  has  been  reduced  to  that  of  a  current-controlled, 
time-invariant,  one-port  device  given  by  [2] 


M(w)=^Ron  +  (l-^)Roff  (3) 


where  w  represents  the  doped  region  of  the  memristor,  D 
the  total  device  length,  and  Ron/Roff  the  lowest  and  highest 
resistance  states  graphically  described  in  Figure  2. 

Equation  (2)  has  been  defined  to  describe  the  velocity  at 
which  w  increases,  meaning  the  device  is  becoming  less 
resistive,  we  use  as  follows  [2] 


dw(t ) 
dt 


Rr 


=  uv^i(t) 


(4) 
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where  uv  is  the  average  ion  mobility  for  the  simplest  case  of 
ohmic  electronic  conduction  and  linear  ionic  drift  in  a  uniform 
field  [2]. 

Utilizing  Ohm’s  law  that  states  that  the  voltage  across  a 
resistor  is  directly  proportional  to  the  resistance  times  the 
current  through  the  conductor;  we  can  obtain  from  equations 
(3)  and  (4)  the  following  relationship  [2] 


The  nonlinear  window  functions  in  Equation  (7a  and  b) 
guarantees  zero  velocity  of  the  doped/undoped  barrier 
interface  described  in  Figure  1  as  w  approaches  either 
boundary,  w=0  or  w=D.  Moreover,  the  differences  between 
the  models  with  linear  and  nonlinear  drift  disappear  when  p 
increases  [5].  The  incorporation  of  the  window  function 
described  by  Equation  (7)  requires  the  redefining  of  Equation 
(4)  as  follows  [6] 


w(t)  =  uv^q(t)  (5) 

and  by  inserting  Equation  (5)  into  Equation  (3),  we  can  obtain 
the  memristance  of  the  system,  which  for  an  Ron  much  less 
than  Roff  can  be  simplified  to  [2] 

M(q)=Roff(l-uv^q(t)^  (6) 


Equation  (6)  describes  the  memristance  of  the  memristor 
system  as  a  function  of  the  charge  q(t). 

Additional  improvements  have  been  proposed  to  the 
aforementioned  memristor  model  to  include  non-linear 
boundary  conditions  [5,6].  The  non-linear  boundary  conditions 
proposed  are  of  the  form 

2  p 

fp(w)  =  1  -  (|  -  stp(-o)  (7a) 

fp(w)  =  1  -  (2^-  1)  P  (7b) 
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Fig.  3.  Memristor  hardware  and  model  fit  correlation  between 
experiment  and  linear  [2],  green  line,  and  nonlinear  memristor  models, 
magenta  [5]  line  and  black  [6]. 


dw(t) 

dt 


=  Wr^i(0/p(w) 


(8) 


IV.  Results  and  Discussion 

Attempts  have  been  made  to  perform  a  model  fit  utilizing 
both  the  linear  [2]  and  nonlinear  [5,6]  memristor  models. 
Figure  3  shows  the  memristor  model-to-hardware  correlation 
fit  utilizing  the  linear  [2],  green  line,  and  the  nonlinear  models 
proposed  by  Joglekar  and  Biolek  [5,6],  magenta  and  black 
lines.  From  the  figure,  it  is  clear  that  the  linear  and  nonlinear 
model  results  do  not  accurately  describe  the  Lissajous  I-V 
characteristic  behavior  of  our  physical  memristor.  Given  the 
lack  of  experiment  and  model  correlation  for  both  linear  and 
nonlinear  models,  we  have  developed  a  simple  empirical 
memristor  compact  model  to  model  and  simulate  the  Lissajous 
I-V  characteristic  behavior. 

One  important  characteristic  of  memristor  devices  is  the  fact 
that  their  present  behavior  is  dependent  on  their  past  state. 
Therefore,  our  simple  compact  model  characteristic  behavior 
is  dependent  on  the  previous  memristor  state  requiring  initial 
conditions.  For  example,  Rmem(t=T)  represents  the  initial  state 
of  the  device  at  time  T,  and  Ron  and  R0ff  represent  minimum 
and  maximum  states.  Assuming  the  device’s  initial  state  is 
Rmem(t=0)=ROn  and  that  in  time  an  input  potential  voltage 
greater  than  Vth  is  present  across  the  device,  we  can  describe 
the  memristance  behavior  follows: 


Rmem  (0  — 
Rmem(t  —  At) 


A tKhleKfl2^Vin^  T ^ 


f  if  ^mem(0  ^  R on 

yelse  Rmem(f)  ~  Ron 


(9) 


where  At=le-3/f  corresponds  to  the  minimum  integral  time 
step  between  memresistance  observations,  f  is  the  frequency 
of  the  sinusoidal  input  voltage,  Th  corresponds  to  the  threshold 
voltage  required  to  enter  the  nonlinear  region  from  the  off 
region,  Khi  and  Kh2  correspond  to  fitting  parameters  used  to 
capture  the  nonlinear  effects. 

On  the  other  hand,  if  an  input  potential  voltage  lower  than 
Vti  is  present  across  the  device,  we  can  describe  the  behavior 
of  the  memristor  device  a  follows: 


Rmem  (0  — 

Rmemit  ~  At)  -  AtKneK^V^)~  W 


f  if  Rmem(.i)  ^  Roff 
[else  Rmeufi)  ~  Roff 


(10) 
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where  Ti  corresponds  to  the  threshold  voltage  required  to  enter 
the  nonlinear  region  from  the  on  region,  and  Kn  and  K12  are 
the  fitting  parameters  to  capture  the  nonlinear  effects. 

Otherwise,  the  state  of  the  memristor  device  remains 
unchanged,  and  the  present  state  of  the  memristor  equals  that 
of  its  previous  resistance  state,  Rmem(t)  =  Rmem(t-At). 

Figure  4(a)  shows  the  correlation  between  the  device 
electrical  characterization  data  and  our  proposed  memristor 
model.  From  the  Figure  4(a),  we  can  observe  that  the  compact 
model  is  able  to  describe  all  regions  of  electrical  operation  of 
the  memristor  device  including  the  high  and  low  conductivity 
regions  and  the  nonlinear  regions.  The  figure  also  includes  the 
linear  model  fit  for  reference  The  compact  model  fitting 
parameters  used  to  achieve  the  compact  model  fit  displayed  in 
Figure  5  are:  Ron=160,  Roff=1200,  Th=0.2,  Ti=-0.35,  Khi=5e6, 
Kh2=-20,  Ku=4e6,  Ki2=20,  an  input  sinusoidal  voltage  of  0.5  V 
amplitude  and  100  Hz  frequency  and  the  initial  condition  that 
Rmem(t)  =  R0ff-  In  addition,  Figure  4(b)  displays  the  time 
dependent  compact  model  results  versus  the  experimental 
results  for  the  memristor  device  hardware  showing  a  good 
model-to-hardware  correlation  in  the  time  domain. 

V.  Conclusion 

We  have  proposed  a  simple  compact  model  that  describes 
the  behavior  of  chalcogenide  based  memristor  devices.  Several 
models  have  been  proposed  in  the  literature  on  the  linear  and 
nonlinear  electronic  characteristic  behavior  of  memristor 
devices;  however,  these  reported  models  do  not  accurately 
describe  the  behavior  of  our  chalcogenide  memristors.  Thus, 
we  have  developed  an  empirical  simple  compact  model  that 
describes  the  electrical  behavior  of  memristor  devices  under 
sinusoidal  input 


Time  (mSec) 

Fig.  4.  Memristor  hardware,  red  circles,  versus  model  correlation  for  the 
linear,  green  line,  and  compact,  blue  line,  memristor  models  (a)  and 
memristor  hardware,  red  circles,  fit  as  function  of  time  versus  and 
compact  model,  blue  line,  (b). 
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